Classifying Vortices in 5* = 3 Bose-Einstein Condensates 
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Motivated by the recent realization of a *^Cr Bose-Einstein condensate, we consider the phase 
diagram of a general spin-three condensate as a function of its scattering lengths. We classify each 
phase according to its reciprocal spinor, using a method developed in a previous work. We show 
that such a classification can be naturally extended to describe the vortices for a spinor condensate 
by using the topological theory of defects. To illustrate, we systematically describe the types of 
vortex excitations for each phase of the spin-three condensate. 
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Over the past few years, there has been remarkable ex- 
perimental progress in the study of spinor Bose-Einstein 
condensates [HM B S, 01 ■ One of the most intriguing ex- 
amples is the recent realization of a ^^Cr [3] condensate 
which has since received substantial theoretical attention 
@j 0, B 0]- Atomic chromium has six electrons in its 
outermost shell which combine into a state with maxi- 
mal spin 5 = 3 according to Hund's rules. These elec- 
trons also give rise to a magnetic dipole moment which is 
much larger than that of the alkalis where dipolar effects 
are typically unimportant. In fact, such dipolar effects 
have been observed in the expansion profiles of ^^Cr [lo| . 
Dipolar effects are important for magnetic traps where 
the spin degrees are frozen out. However, when purely 
optical traps are used it has been argued that the spin 
exchange interaction will overwhelm the dipolar interac- 
tion . Due to such spin interactions, a dilute gas of spin 
three bosons will interact according to the pair potential 

Vp{ri - ra) = <5(ri - r2)(ffo^o + 92V2 + giV2 + 56^6 )(1) 

where Vs projects into the state with total spin S and 
gs = 4Trh?as /m where as is the scattering length cor- 
responding to spin S. For the spin-three problem, this 
potential gives rise to a phase diagram (as a function 
of the possible values of the scattering lengths) of great 
richness as was revealed in [1, 0] • In addition to having 
such rich phase diagrams, spinor condenstates in gen- 
eral and the ^^Cr condensate in particular will have in- 
triguing topological excitations. Topological excitations 
for spin-one, spin three-halves, and spin-two conden- 
states have received considerable attention in the liter- 
atu re [ll , l3, 13, li , 3 lB| since the pioneering work of 
Ho [iVj and Ohmi and Machida Also, recently topo- 
logical spin textures have been observed in '^^Rb atoms 
in the spin-one hyperfine state [1]. 

In a previous publication [l^ , we have developed a way 
of classifying spinor BEG and Mott insulating states and 
illustrated the method with an application to the spin- 
two boson problem. More specifically, we showed that it 
is natural to identify a spin S state 1-0) = Sf=-s^al'^) 
with 2S points on the unit sphere {9,4>). In fact, there 
is a one-to-one correspondence between a spinor and a 



"reciprocal spinor" up to phase and normalization 



(2) 



where C = tan(6'/2)e** is the stereographic image on 
the complex plane of the point [6, 0) on the unit 
sphere. The set of such complex numbers {Q} are 
explicitly constructed by finding the coherent states 

IC) = Ealo \/OC"|5' - a) which are orthogonal to the 
spinor states 10). The characteristic equation f.^,{C) = 
{ip\C) will then be a polynomial in ( of order 25* and its 
25' roots projected on the sphere will immediately give 
the symmetries of the spinor. Spin rotations which leave 
this set of points invariant will take a spinor Aa into 
itself modulo an overall phase which can be computed 
through geometrical considerations. Thus this method 
gives a complete description of the spinor state. 

In this work, we show that such a method is also use- 
ful for classifying the spinor vortices. To illustrate, we 
focus on the spin-three problem, motivated by the re- 
cent experimental development [J|. We first reproduce 
the phase diagrams for spin-three condensates obtained 
previously in [0, 0] . Then we find the symmetries of each 
phase (see Fig. [T|) by applying the method developed in 
[19| . We then systematically describe the types of topo- 
logical excitations for such a system using the topologi- 
cal theory of defects [2^. Such a scheme enumerates all 
possible types of vortex excitations, though it does not 
determine which are energetically stable against dividing 
into multiple vortices. Finally, we will briefly discuss the 
resulting microscopic wave functions for the vortex exci- 
tations. Such understanding of vortices will be relevant 
for rotating ^^Cr condensates in optical traps. 

In general, the topological vortex excitations (line de- 
fects) of an ordered medium are given by the first ho- 
motopy group 7ri(i?) (the group of closed loops) of the 
order parameter space R. More specifically, the conju- 
gacy classes of 7ri(i?) correespond to the vortex types. 
Moreover, the way such vortices combine is determined 
by the multiplication table of the conjugacy classes. For 
a pedagogical review of the topological theory of defects, 
see 120| (whose notation we will adopt) . A theorem from 
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algebraic topology (see e.g. [2l[) which is used as a reli- 
able tool for computing the first homotopy group is the 
following. Let G be a simply connected group and H be 
a subgroup of G. Then 

7ri(G/F) = H/Ho (3) 

where Hq are the elements of H which arc connected 
to the identity by continuous transformations. To ap- 
ply this theorem to our system, we need to take G to 
be a group of operators that act transitively on the spin 
groundstates and H to be the subgroup of G which leave 
the particular spinor state we are considering invariant. 
Then R has the same topology as G/H, so the theo- 
rem applies. The first guess would be to use elements 
like e*"e~'^ "'' (where S is the spin-three matrix and n 
is a unit vector) which act directly on the spinor states 
and superfluid phase. However, these elements form a 
representation of 5*0(3) (g) C/(l), which is not simply con- 
nected, so Eq. ^ will not apply. Thus we need instead to 
take the simply connected group G = SU{2) (g) T where 
T is the group of translations in Id with elements sat- 
isfying the property T{x)T{y) = T{x + y) (which will 
be simply connected). This is the covering group of 
5*0(3) (X) C/(l). In this work we will use the represen- 
tation of SU{2) by Pauli matrices, thus we describe its 
elements in the form e"'^^®. We then use the homo- 
morphism ip : SU{2) ig) T 50(3) (g) U{1) given by 

ip (T(a;)e-'T^®) = e^e^'^ "® which "lifts" our group 
acting directly on the spinor states into the larger simply 
connected group, G. The subgroup H C G will be the 
collection of elements h G G which leave the particular 
spinor under consideration invariant. Once H is con- 
structed by such a method, Eq. ([3]) can then be directly 
applied to construct the first homotopy group. 

As we will see, the majority of the phases in the spin- 
three phase diagram have only discrete symmetries. For 
such states, Hq is the trivial group containing one el- 
ement. For this situation, our computational theorem 
reduces to 

TTi{G/H)=H. (4) 

Therefore, constructing tti for such states becomes 
straightforward: One first finds the rotations in 50(3) 
which leave invariant the points on the unit sphere {Ca} 
given by our classification scheme [l^ . For this purpose, 
it is helpful to note the possible finite subgroups of 50(3). 
They are the cyclic groups 0„, the dihedral groups £>„, 
the tetrahedral group T, the octahedral group O, and 
the icosahedral group Y (see Table H]). One then lifts this 
subgroup of 50(3) into the twice-as-large group SU{2). 
On top of this "spinor scaffolding" one can impose a su- 
perfluid winding x, but the allowed values of x are re- 
lated to the spin rotation chosen. The winding is given 
in the form x = —X+2T:k where the integer k is arbitrary. 



The offset A is fixed and can be determined geometrically 
through the relation 

A = (r - 5)9 (5) 

where O is the angle of rotation, r is the multiplicity 
of the spin-root on the axis of rotation, and S is the 
total spin (see appendix for proof). This information 
on the superfluid phase determines the translation group 
elements T(x). This provides a direct way of constructing 
H and therefore classifying the possible vortices which we 
illustrate below for the spin-three problem. 

Now we move on to discuss the possible phases of the 
spin-three condensate in the absence of an external mag- 
netic field. The interaction hamiltonian determined by 
the pair potential Eq. ([1]) is given by 

Hint = -;^'^gsa^aa\{al3\Sm){Sm\a' (3')aa'ai3i (6) 

S.m 

where ajj creates a zero-momentum boson in the spin 
a state, (a/3|5m) are Clebsch-Gordan coefficients, and 
the sums over the Greek indices are implicit. Note that 
this interaction energy alone is sufficient to determine 
the spinor ground state since the rest of the hamiltonian 
will not have a spinor preference. Using the variational 
state It/i) = (Acajj) |0), where the coefficients are 
normalized A*^Aa = 1 the interaction energy to leading 
order in N can be found by replacing Oq, by the scalar 
^/N Aa giving the interaction energy of 

E,,jN' = ^Y.9sPiS) (7) 

s 

where 

P{S) =Y,\A:A;{aP\Sm)\\ (8) 

rn 

We numerically minimize this energy as a function of 
the seven complex variational parameters Aa under the 
constraint of normalization A^A^ = 1 for all possible 
values of the scattering lengths, and reproduce the phase 
diagram given in 0, 0| ■ We give explicit forms of the wave 
functions up to 5*0(3) rotation which are summarized in 
Table U where, for comparison, we use the same notation 
as was used in A subtle difference is that we find 
the additional state HH which is degenerate with the 
phase H reported in Q. It is important to note that 
the unspecified angles 77 and x will take unique values for 
particular regions in the phase diagram. That is, there 
is no mean-field degeneracy as was reported for the spin- 
two case figj . 

We now apply the classification scheme developed in 
[igt to give the symmetries of the possible spinor con- 
densates, and the results are given in Fig. [TJ Shown are 
reciprocal spin-three spinors which arc six points on the 
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TABLE I: Wave functions \tp) = (^Qa^) |0) for the possible spinor phases corresponding to global minima of ^ for 
particular scattering lengths up to SO{3) invariance. Also shown is the subgroup of SO(3) corresponding to the symmetries of 
the wave function. The last column gives the number of conjugacy classes of tti for each spin state, which gives the number of 
possible vortex excitations. The "xZ„" in this column reflects the fact that an arbitrary multiple of 2n can always be added 
to the winding in the superfluid phase. 
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FIG. 1: Reciprocal spinors. Points on the unit sphere (gray 
spheres) having the transformation properties of the wave 
functions given in Table |I] The small white spheres are put 
at the origin as a visual aid. Points which correspond to de- 
generate roots are marked. 



unit sphere which make the symmetries explicit. All of 
the phases except for F and FF have no continuous sym- 
metries. As was explained in [1, 0], all of the scattering 
lengths for ^^Cr arc known except qq. This gives a line in 
parameter space of possible ground states for such a con- 
densate where the phases A, B, and C are the candidate 



ground states of the ^^Cr condensate. 

We will now discuss the topological excitations of each 
phase. For each case we will take G = T ® SU (2) which 
will act on the spinor states via the homomorphism ip. 
The goal is to find the isotropy group H which leaves the 
particular spinor state under consideration invariant. We 
start with the two states having continuous symmetries, 
namely F and FF. For state FF, taking the state to be 
oriented along the z-axis, we find that the only operators 
leaving the spin state invariant are ^i-SQ^-^S^Q with S = 
3. Note that the phase factor e**^® is necessary to cancel 
the phase from rotation. Immediately, we find that this 
gives the isotropy group 



H = {r(5e + 27rn)e"**® : e e 

Ho = {r(S'e)e"*'^® : e e M, } . 



,nez| (9) 
(10) 



The first homotopy group is then computed 
7ri(FF) = H/Hq = lis consistent with that in 



to be 

Carrying through the similar analysis for state F, we 



find 7ri(F) = liiis-i)- Note that since these groups are 
abelian, the conjugacy classes are the original group ele- 
ments giving. The number of conjugacy classes for these 
cases are 25" = 6 (FF) and 2(5*- 1) = 4 (F). The number 
of possible types of vortex excitations is smaller by one 
here since the identity corresponds to the absence of a 
vortex. As noted before, this is quite a nontrivial result. 
This is in contrast to the spinless case where there are in- 
finitely many types of vortices corresponding to different 
windings in the superfluid phase as the vortex is circled 
(with TTi = Z). The underlying reason for this is that ro- 
tating such a spinor along its symmetry axis is identical 
to changing the superfluid phase. 

We now move on to consider the remaining phases 
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which have only discrete symmetries. For these situa- 
tions, as mentioned before, Hq is the trivial group, so we 
can apply Eq. (g]) directly. We will start with the phases 
having the smallest isotropy groups, namely phases C 
and G. For phase G, with the spinor wave function writ- 
ten as in Table HI the only symmetry is a rotation by tt 
about the z-axis. Such an operation will not change the 
phase of the spinor, i.e. e~^^'''^\ip) = \tp). We then find 
for the isotropy group 

7ri(G) ^ H = {±To,±icr^} (g) {T{2TTn) : n e Z} . (11) 

where tq is the 2x2 identity matrix and we note that 
^±i^ir _ ±ifj^_ This group is abelian, so the conjugacy 
classes are just the group elements themselves. Note that 
the translational part of H merely tells us that an arbi- 
trary winding in the superfluid phase can be imposed on 
any vortex state. In the following, we will suppress the 
translational part of H, keeping this in mind. The phase 
C will have a similar isotropy group, but we note here 
that rotation by tt about the symmetry axis will give a 
phase of = — 1. 

Phases E, H, and HH have n-fold rotational symme- 
tries about a particular axis. Phases E and HH have 
a three-fold rotational symmetry about their symmetry 
axes, and the first homotopy group is found to be 



7ri(E,HH 



2 3 : TT, 



1 



,2,3}. (12) 



On the other hand, phase H has a five-fold rotational 
symmetry, and its first homotopy group is similarly 



7ri(H)=={ 



±e- 



(13) 



As before, these groups are abelian. Therefore, the types 
of possible vortices will correspond to rotating the poly- 
hedra about the axes of symmetries n times as the vortex 
is circled. In addition, the phase winds by — A+27rfc where 
k is an arbitrary integer and A = 0, — for the 

phases E,H, and HH respectively. 

Phase B is the first phase we encounter that has a 
nonabelian isotropy group. This phase has a rotational 
symmetry of tt about the z-axis. In addition, it also has 
the symmetry of a rotation about two axes in the xy 
plane, namely {nx,ny) — -^(1,±1). The isotropy group 
for this is 



^i(B) 



±To,±iaz,±—i={a-x ± (Ty) 
v2 



(14) 



Since this group is nonabelian, the conjugacy classes are 
not just the original group. This group is isomorphic to 
the quaternion group having eight elements. The conju- 
gacy classes corresponding to each possible type of vortex 
excitation are: 

Co = {to} , Co = {-To} , = {±ia,} , (15) 



Cxy - <j ±-J=i<Tx + (Ty) 



Cxy = { ±-J={(yx - (Ty) 



Phases A and D have larger isotropy groups; we 
will not write out explicit representations here, but will 
merely state the results. For phase A which has the sym- 
metry of the hexagon as seen in Fig. [1] there is a six- fold 
rotational symmetry about the z-axis. In addition to this 
there are six axes in the xy plane for which there will be 
a two- fold rotational symmetry. This will make 7ri(A) 
have 24 elements, which will have 9 conjugacy classes. 
Phase D has the symmetry of the octahedron. For this, 
7ri(C) will have 48 elements and 8 conjugacy classes. All 
of these results are summarized in Table H] 

Now that the classification scheme for vortices based 
on homotopy theory has been completed, we will briefly 
consider the resulting microscopic wave functions. We 
can express the order parameter explicitly in terms of the 
spinor as il^ai^) ~ ilj{Y)Aa{v) and write V'(r) = |?/'o|e'^(''\ 
Note that the amplitude of '0(r) is taken to have con- 
stant magnitude which is realistic provided we are suf- 
ficiently far from the vortex core. For Aa{v) it will 
prove useful to use the spinor-ket notation and write 
\A{r)) = e~*^ "'^('')|Ao) where \Ao) is the reference state 
which can be taken as one of the spinors given in Ta- 
ble [D In this equation, n is the axis about which the 
spinor is rotated as the vortex is circled and 0(r),/3(r) 
give the spatial dependence of the superfluid phase and 
spinor rotation around the vortex. These functions must 
satisfy the requirement that the original wave function 
must be recovered after the vortex is completely circled. 
For instance, suppose the spinor state under considera- 
tion is invarient under a rotation by angle 0. Then a 
vortex at the origin corresponding to this symmetry is 
given by 0{y) = (— ^ + k)ip and /3(r) = where ip is 
the azimuthal angle, k is the winding number, and A is 
given by Eq. 
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Note added: After this work was completed, we be- 
came aware of a work by S.-K. Yip (2^ which considers 
a similar problem of vortices in spinor BECs, but uses a 
different approach. 



APPENDIX: PROOF OF EQ. (HJ 

The geometrical reciprocal spinors determine the rota- 
tional symmetries of the spinors. Associated with each 
such symmetry is a phase factor which we will calculate 
here. Suppose the rotation through an angle Q about 
the axis n is a symmetry of the reciprocal spinor of 
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Then 



(16) 



where A is given by Eq. 

To prove this, it is convenient to choose (without loss 
of generahty) the axis of rotation to be parallel to the 
z-axis: n = z. Then by comparing the components of 
the spinors on the left- and righthand side of Eq. ([16]) we 
see that 



-iaQ 



(17) 



This equation implies that the non-zero components of 
the spinor must satisfy 



A; 



-aQ (mod 27r) 



(18) 



Let us compare this result to the multiplicity of the 
root at the north pole, which corresponds to C = 0. The 
characteristic equation is 



MO 




.C = 0. 



(19) 



Q=0 



If the root at zero has multiplicity r, then the coefficients 
of the first r terms ■ • ■ , ^re all equal to zero, and 
the r-\- 1*'' coefficient is nonzero. This corresponds to the 
r -|- 1*^' component of the spinor being nonzero: 



As-r ^ 0. 



(20) 



Substituting m = S — r into the condition for the nonzero 
components Eq. gives the result Eq. ([5]) and the 

proof is complete. Note that our discussion has been for 
a rotational axis aligned with the z-axis, but the result 
Eq. ([5]) has a conceptual formulation which does not de- 
pend on the coordinate system. 
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